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Abstract 

We compute how the initial energy density and produced gluon, quark and antiquark 
numbers scale with atomic number and beam energy in ultrarelativistic heavy ion 
collisions. The computation is based on the argument that the effect of all momentum 
scales can be estimated by performing the computation at one transverse momentum 
scale, the saturation momentum. The initial numbers are converted to final ones by 
assuming kinetic thermalisation and adiabatic expansion. The main emphasis of the 
study is at LHC and RHIC energies but it is observed that even at SPS energies this 
approach leads to results which are not unreasonable: what is usually described as a 
completely soft nonperturbative process can also be described in terms of gluons and 
quarks. The key element is the use of the saturation scale. 
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1 Introduction 



In ultrarelativistic heavy ion collisions, the number of produced gluons and quarks 
with pt larger than some cut-off po, Naa{po, V^), increases when po decreases, 

when the size of the nucleus increases (~ da/aind ~ ^^/-^a ~ A^^^) and when ^/s 
increases (since the small- a; enhancement of the distribution functions becomes 
more effective). Shadowing will decrease the number Q], but NLO corrections will 
increase it. At sufficiently large cut-off po ^ Aqcd, the system of produced gluons 
and quarks is dilute and usual perturbation theory can be applied. However, at some 
transverse momentum scale po = psat the gluon and quark phase space density saturates 
and one does not expect further increase. In this case one may conjecture that 
the amount of produced at Psat gives a good estimate of the total Et produced in 
an average nucleus-nucleus collision: gluons with px ^ Psat carry lots of E^ but are 
very rare, whereas gluons with pt ^ Psat are numerous but carry little Et- Doing 
the computation at psat gives an estimate of the effect from all scales, both above and 
below Psat- 

In this note, we shall work out the quantitative consequences of the above conjecture 
for the amount of Et initially produced in an A + A collision at some large ^/s. The 
key element is the definition and computation of the saturation scale. For this, one 
first computes the number Naa{po, y/s) of produced quanta including perturbatively all 
gluons and quarks with pt > Po- A saturation criterion is formulated and the solution 
of this gives psat- Evaluating NAA^Psat, a/s) and E^'^i^psat, a/s) then gives the initial 
values. All quantities can well be fitted to a scaling law of the type {y/s is in units of 
GeV in the formulae) 

CA'^i^sf, (1) 

possibly with an additional term constant in s. 

This computation at ps^t only gives initial values at the proper time = 1/psat and 
not yet experimentally measurable quantities. To get these, one has to trace the system 
through the entire set of expansion, hadronisation and decoupling stages. We shall also 
discuss this assuming kinetic thermalisation of the gluonic component and adiabatic 
expansion. 

One is accustomed to considering the average events at SPS energies (yi< 20 GeV) 
as essentially soft, stringlike. It is somewhat surprising that this entirely hard, per- 
turbative, treatment will be seen to give numbers which are not unreasonable even at 
these relatively low energies. This phenomenon has also been pointed out in |Q TU 



m 



the framework of a parton cascade model. From the point of view of our study the 
crucial factor is the use of the saturation scale as the scale of dominant processes. 

There are clearly many refinements (or objections) one can suggest to this approach. 
The most obvious is that one should try to spread out the relevant momentum scales 
around Psat or, in time, around Tj = 1/psat- Further, there certainly will be some 
entropy production during the expansion stage - that due to chemical equilibration 
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has been estimated in |Tl|. We will discuss these briefly, but the key issue anyway is 
the use of Psat as the dominant scale. 

If Psat is very large, well into the perturbative regime, there are transverse scales 
Pt ^ Psat which are perturbative but for which transverse phase space occupation 
numbers are ^ 1. A classical field description is then appropriate |12|-[^, in analogy 
with the computation of the rate of baryon number violating reactions in electroweak 
theory |pT|. 



2 Initial values 

The computation is technically a standard pQCD one and uses the formulas in Sec- 



tion 2 of |2^. The parton distribution functions are from |2^ and nuclear effects to 
them are implemented using the EKS98-parametrization @]. The corresponding NLO 
corrections to Et within the appropriate acceptance region {\y\ < 0.5, all azimuthal 
angles around the beam direction, and px > Po) have been computed in p, see also P]. 
This is an infrared safe computation. For, say, po = 2 GeV and choosing all the scales 
as pt, one finds K = 1.7 &t ^/s = 5500 GeV, K = 2.3 at ^/s = 200 GeV and K = 1.9 at 
^/s = 20 GeV. In principle, the K-factor should be computed at the saturation scale, 
but in view of this small range of variation we shall take K = 2. Anyway the NLO 
corrections to N cannot be computed in an IR safe manner without introducing some 
additional regulator. 

The calculation now proceeds as follows. First compute the average number Nk{po), 
k = g,q, q, of quanta produced in an average central A + A collision within \y\ < 0.5 
(for brevity, the fact that various quantities may depend on A, centrality, i/i or pq is 
not explicitly marked in what follows). The saturation limit is then obtained from the 
equation 

iV(po)= E NM=pIR% (2) 
where Ra = 1.12^^/3 rpj^jg 

expresses the fact that at saturation A^(psat) quanta 
each with transverse area /pls^t the whole nuclear transverse area 'kR\. No nu- 
merical or group theory factors nor powers of Qs are included; these are anyway 0(1) 
unless one discusses a parametric weak coupling limit — 0. All parton flavours are 
included, though at Pq = Psat gluons clearly dominate even at lowest energies. The way 
Eq. (0) determines p^at for A = 208 and ^/s = 17, 200, 5500 GeV is shown in Fig. |l|. 
The numerical values computed for various A and -^i form a family of curves shown 
in Fig. ^ (normalised to give final charged multiplicity). They can be well fitted by 
(the points with ^/s > 200 GeV are used in determining the parameters, though the 
fit agrees well also with the points computed at lower y/s) 

iV. = Ar(po = Psat) = 1.383 A°-922(^)0.383^ (3) 
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Psat = 0.208 GeVA°-i2«(v/^)°-i^\ (4) 
Ps.,Ni = 0.288 GeVAi o^O(yi)°-5^^ (5) 

With given assumptions {K = 2, certain set of parton distributions, certain method 
of implementing shadowing, certain scale choices, etc) we require the accuracy of the fit 
to be within 1% at the fitted region and quote the parameters with corresponding ac- 
curacy. Clearly there is a larger error related to the underlying theoretical assumptions 
(e.g. the ii'-factor) and to the whole starting point: the use of Psat as the dominant 
scale. 




Po [GeV] 

Figure 1: The average number of initially produced QCD-quanta with pT > Po and \y\ < 0.5 
as a function of the lower limit pq for central Pb-Pb collisions at y/s = 5500 (LHC), 200 
(RHIC) and 17 GeV (SPS). The saturation scale psat for A = 208 is given by the points of 
intersection of the dashed curve "saturation" (p^R^) with the curves N{pq). 

Secondly, with the psat so obtained, the initial values ETi,k and their sum can be 
computed. The results for the sum are plotted in Fig. ^ and they behave as 

En = 0.386 GeV Ai-°^=^(v^)°-^^^ (6) 
= 1.34A~°-°°^(v/^)°-°>.,,iV,. (7) 

By comparing (|^) and (|^ one observes that approximately Eti ~ const Psat-^- Ac- 
tually the dependence on ^/s and, in particular, on A is so weak that one may question 
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whether the deviation is a genuine physical effect. To interpret the result, assume first 
that precisely 

ETi = Cip^^tNi, (8) 

with Ci = constant. Then, using the initial volume estimate Vi = i^R\/psa.t and the 
saturation relation Ni = Ps^.t^A o^i^ obtains 



Psat, (9) 



Ni 1 . 

= T7 = -Psat, (10) 
Vi IT 

- = CiPsat. (11) 

rii 

Since the gluons are by far the dominant component, even at the SPS, let us convert 
ej to a temperature including only them, i.e., writing e = IGvr^/SO ■ T^. Then 

Comparing (|lTD and (0), the boson gas thermal relation e/n = 2.70T (=3.15T for 
fermions) is seen to hold if Cf = 2.7^30/(167r^), Ci = 1.48; very close to what the 
numerical calculation gives for it (see Eq. (0)). 




Figure 2: (a) The initial Et in |y| < 0.5 in a central A + A collision with A = 
12, 32, 64, 136, 208 as a function of ^/s. (b) Decomposition of the initial Et into gluon, quark 
and antiquark components for A = 208. 
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Quantity 


C 


a 


h 


SPS 


RHIC 


LHC 






Q29 


383 


598 


1440 

X X xw 


5140 


n JGeV 

i-'sat/ » 




128 


IQl 


73 


1 13 

X . X 


2.13 


ETi/GeV 
e,/(GeV/fm3) 
Hi ■ fm^ 


0.386 


1.043 


0.595 


600 


2360 


17000 


0.103 


0.504 


0.786 


16.1 


98.2 


1330 


0.370 


0.383 


0.574 


16.0 


59.8 


401 


Ti/GeV 


0.111 


0.126 


0.197 


0.39 


0.62 


1.19 


Etc 


3.477rc 


0.917 


0.398 


276 


692 


2600 



Table 1: Summary of the values of the scahng parameters C,a,b in quantity 
together with numerical values for A = 208 at SPS, RHIC and LHC (y/s = 
GeV). The numbers for Etc (Eq. (H) are for = 0.18 GeV. 



20, 200, 5500 



Including the small deviations from 
En 



Hi 

T- 
rii 



CiPssitNi one finds, similarly: 

-- 0.103 GeVfm-=^A0-50^(v^)0-^8^ 
0.370 fm-3A0-383(y^)0.574^ 



0.111 GeV A°-i26(y^)0.i97 
2 52^-o.oo5^y^)0.oi5y^_ 



0.534A 



-0.002 



(v^) 



0.006 



Psat; 



(13) 

(14) 

(15) 
(16) 



One observes that the ideal boson gas thermalisation relation e/n = 2.70T holds al- 
most independent of A and ^/s^. the system satisfies this criterion of thermalisation 
right at production. Note that this applies to some degree even at the SPS. The scal- 
ing parameter values together with some relevant numerical values are summarised in 
Table 

It is also illustrative to consider the results for g,q,q separately. The A and y/s 
dependences of the different components are almost the same but to account for the 
growing importance of the valence quark component towards smaller energies, one also 
has to include a constant term in the fit. For simplicity, we first make a four paramenter 
fit A''[C{^/sY+D] to the total initial muhiplicity Ni for y/s > 40 GeV. With the powers 
a and b so fixed, we find C and D separately for each component. The result is 



A°-922[i_o65(v^)°-^°^- 0.028], 
A°-'^^[0.021{V^f^'' + 0.778], 
A°-922[o.o37(v^)°-^°^ + 0.287], 



(17) 



and for the initial transverse energy 

ETg = Ai-°^^[0.341(v^)°-599_ 0.398] 



:i8) 
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= Ai-°^^[0.0137(v^)°-^^9 + 0.494], 
Erg = Ai-°^^[0.0153(v^)°-''' + 0.153]. 

Note that the simple parametrisations above are only meant to reproduce the behaviour 
of the different components, not the small difference Ng — Ng. They cannot be used to 
compute, e.g., the net quark number at the LHC, since (A^^^ — Nq)/Ng < the accuracy 
of the fits. 

3 Expansion stage 

We have now the initial values at = 1/psat- To compare with experiment one has 
to follow the evolution of the system through expansion in the QCD plasma phase, 
phase transition, expansion in the hadronic phase and decoupling. Now that one has 
computed from perturbative QCD that the initial state very nearly satisfies the kinetic 
thermalisation condition e/n = 2.7 T, there is some justification in assuming that also 
further expansion is locally thermal, i.e., entropy conserving. There is a large number 
of dissipative effects one can think of, but adiabatic boost invariant expansion is the 
baseline here. 




r[fm] 

Figure 3: Proper time dependence of energy density during longitudinal expansion for LHC, 
RHIC and SPS with the initial values given in Table |. 

Numerical examples of adiabatic expansion are shown in Fig. ^ for the initial one- 
dimensional stage during which e(r)r'^/^ = constant. We start the expansion at the time 
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1/Psat, although at lower energies this may be smaller than the transit time 2Ra/'^ = 
AlRaI \fs of the nuclei. From the scaling formulas it actually follows that IR^j^ < 
1/Psat for 

> 6.81 A^-^™ (= 143 for A = 208). (19) 

We remind that the whole production process is effectively described as taking place 
at the time 1/psat- 

Later on, at r ~ ^/SRa, the system will enter a regime of longitudinal and transverse 
expansion, which has to be studied numerically |2^. In Fig. ^we see that at the LHC 



the transverse expansion effects are likely to set in already in the plasma phase. 

For the multiplicity the situation is simple: for an ideal system of bosons and fermions 
the total entropy is S" = 3.60Nb + 4:.20Nf and this is constant per unit rapidity, as 
long as the expansion is boost invariant. (By definition, all quantities are always for 
\y\ < 0.5.) Thus, initially. Si ~ 3.6A^j. For the final hadronic gas Si = S'/>4iV/ so 
that Nf = O-dNf. the number of final state hadrons is to up to 10% corrections equal 
to that of the initially produced partons (gluons) at the scale psat- The multiplicity 
prediction Nch = 2/3 * 0.9Ni thus is directly given by Eq. (|^) and is shown in Fig. ^. 

Since NLO corrections to Et can be computed in an IR safe manner, one might 
prefer to compute the final entropy from the initial through the conversions ETi 
ei ^ Ti —>■ Si = Sf. By making use of the fit (p^, we obtain the following scaling law 
for the multiplicity of charged particles per unit rapidity: 

Nch = ^l^ = ll.WA'-^'iV-sr". (20) 



In view of the observed initial near thermalisation (Eq. (|T6|)) this necessarily leads to 
the essentially same scaling in A and ^/s as that obtained from Ni in Eq. (^). The 
comparison is shown in Fig. ^ 

For Et the situation is more complicated. Et per unit rapidity decreases by the 



factor T/Ti due to work done against expansion On the other hand, the 

transverse expansion at later times will compensate for this effect. 

To estimate the final Et, we extract an approximate scaling law between the final 



average transverse momentum (pt) and the final multiplicity Nj from Fig. 7 of Ref . |^ : 



{pT)/GeV = 0.39 + 0.061 \n{Nf/A). (21) 

By substituting our total multiplicity Nf, the final transverse energy can then be 
computed as 

ETf = Nf{pT). (22) 

The transverse energies so obtained are shown in Fig. |^. Using Eq. (^) we get the 
scaling in A and ^/s^. 

Etj = 0.46A°-92^yi)0.4o^^ _Q_Q^21nA + 0.0611nyi]. (23) 
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Figure 4: The number of charged particles per unit rapidity for A = 12, 32, 64, 136, 208 as a 
function of ^/s, computed from Eq. ( |20| ) (sohd hues) and compared with 2/3 * 0.9A'"j (dotted 
hues). 

On the other hand, a simple estimate of the transverse energy E^c at the end of the 
plasma phase is obtained by neglecting the transverse expansion during the plasma 
phase. This gives 

Etc = {TJT,)Et. = 3.477T,A0-9i^(v^)0-398 (24) 
= 2.51T,A-o-°°5(v/^)0-°i5iV,. 

Using here = 0.18 GeV, E^c approximates the final Ejn computed above very well. 
The comparison is shown in Fig. 

Eq. (p^) expresses a trivial fact: as long as the system is thermal, the energy per 
particle is ~ 2.7T - possibly even larger if also quarks enter chemical equilibrium. 
Evolution in the hadronic phase until decoupling will lead to a further reduction by 
the ratio TdecoupUng/^c^l, but this is compensated for by the development of flow. Thus 
the net effect of the hadronic phase can be expected to be small, and Etj ~ Etc- 

In Fig. we show the difference between the initially produced Etx and the final, 
measurable, Exf for A = 208. At the LHC the transverse energy drops by a factor ~ 6 
and at RHIC by a factor ~ 3.5 due to pdV work done against expansion. 



8 




Figure 5: (a) The average final Et computed from Eq. (^) (solid lines) and from Eq. (|2^ ) 
with Tc = 0.18 GeV (dashed lines) for A = 12,32,64,136,208 as a function of ^/s. (b) 
Comparison of the initial and the final Et for A = 208 as a function of ^/s. 

4 Discussion 

In the saturated scaling limit (ex ~ p~t) one expects the following A-scaling rules: 
First from Eq. (g), N ~ v4'^/VPsat = (Psati?yl)^ one obtains Psat ~ A^^^ and N ~ 
A. Further, Eti ~ Psat^ ~ A"^^^. Due to shadowing the observed A-exponents are 
somewhat less than these numbers, for Ni 0.922 instead of 1 and for Eti 1.043 instead 
of 1.167. Anyway it is quite interesting to note these rather small exponents for this 
hard process. Note that for a fixed momentum scale without saturation the growth 
with A is much more rapid: iVj ~ E^i ~ A^^^. 

The growth with beam energy is rather rapid, the power of y/s is 6 ~ 0.19 for psat, 
b ^ 0.38 for Ni and about the same for the final transverse energy, depending on the 
flow effects. A constant valence quark component can also be identifled. The power 
of ^/s depends on the phase space available and on the small-x increase of the gluon 
distribution function. 

The decrease of transverse energy in adiabatic expansion (Fig. ^) is rather large, 
a factor ~ 6 (3.5) at the LHC (at RHIC). Provided that one can control flow effects 
numerically, a measurement of Exf will give further constraints on the degree of ther- 
malisation in the system. 

An interesting property of saturation calculations appears when considering them at 
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SPS energies. Usually a heavy ion collision at the SPS, y/s<20 GeV, is regarded as 
entirely soft, nonperturbative, dominated by stringlike effects. Here we have gone to the 
other extreme and have treated it in terms of the gluon and quark degrees of freedom 
at a very specific scale, pt = Psnt- This scale is only 0.73 GeV for Pb+Pb collisions at 
the SPS, and perturbation theory at this scale, although a;s(psat)/7r is clearly < 1, is 
not expected to be very accurate (though it works well for, say, r hadronic decays with 
rrir = 1.8 GeV ^6|). However, the numbers one obtains (Table j^) are not very far from 
the measured ones. The charged multiplicity per unit y (actually rj), is experimentally 



2^ about 400 for Pb+Pb at = 17 GeV while we obtain N^h = 330.. .340. For S+S 
at ^/s = 20 GeV, the experimental value is N^h ~ 50 while we get Net ~ 62. ..64, 
again not dramatically different. The average for Pb+Pb is observed to be about 
400 GeV [^], somewhat larger than the numbers in Table |T]. We have not in any way 
attempted to improve the quantitative agreement with present data. Maybe some kind 
of duality is at work here. 

It is also notable that even at the SPS the system at psat is dominantly gluonic: at 
y/s = 17 GeV the gluons take 65% of ETi and 73% of Ni. A related fact is that even 
at the SPS the thermalisation condition e = 2.7Tn is initially nearly satisfied. 

The RHIC numbers {A = 208) are ^1300 for the total muhiplicity dN/dr] (870 for 
A^ch from (pOD) and pa660 GeV for dEx/drj and will soon be measured. However, one 
can compare them with some historical numbers. In [|^ the values psat = 0.95 GeV and 
1300 estimated for the multiplicity are very close to the ones computed here, but the 
discussion of Et included no p dV work. 

For the LHC, we get the charged particle multiplicity dNch/dr] ^ 3300 and dET/drj ^ 
2900 GeV. These are some 25% larger than the LHC-numbers in Ref. p2|]. Now 
different, largely compensating, effects are included: shadowing and a slightly larger 
value of Psat decrease the numbers while K = 2 instead of K = 1 increases them. 

Finally, some remarks are listed: 

• The saturation criterion (^ has been presented here in a very simple form. Fun- 
damentally, it is the result of a complicated dynamical computation containing 
numerical and group theory factors and powers of Qs The precise value of 
all these will affect the final numbers. Taking the net effect to be = 1 in (0) is 
a geometric estimate based on the uncertainty principle and there is no way of 
giving a controlled estimate of the error. 

• The saturation criterion could also be formulated by applying it to the pt- 
distribution, calculating it perturbatively for pt > Psat and assuming that it 
is, say, constant for pt < Psat- This would lead to somewhat different numbers, 
but the error is like that under the previous point. 



To relate the present study to Ref. [30|, it is perhaps interesting to compute the 



screening mass mg = gT at T = Tj. Using Eq. ( [T3| ) we get mg = gs{Psa.t)Ti 
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1.6 In ^''^(psat/AqcD)^ °'°°^(v^)° °°^Psat- With the saturation scales in the region 
Psat = 1---2 GeV, we get nig ^ Ps^t, since p^at ~ O.STj. 

The saturation criterion is formulated for quanta produced in the collision. The 
wave functions of the initial nuclei need not be saturated for the values of A and 
^/s considered here. 

The starting point of these calculations of initial production of QCD quanta in 
ultrarelativistic central heavy ion collisions was that it is enough to consider one 
momentum scale, pt = Psat- The magnitude of the various quantities was then 
computed using methods appropriate for px > Psat- One can, as well, attempt 
to use the range pt < Psat- This is a region of QCD in which genuinely new 
dynamical methods have to be developed. Numerical calculations also there are 
under way |TB|, However, whatever the method applied, it anyway applies 



only to the initial production (1/psat < t < 1 fm), leaving the treatment of the 
entire further evolution of the system open. 

There are many interesting and relevant questions concerning the subsequent 
flow, after initial production. The key issue is thermalisation, that of quarks, 
in particular |TT|, One could formulate the problem with all degrees of free- 
dom and with some suppression factors (fugacities |^) to account for the initial 
deficit of the quarks. To avoid this complication but to partly account for the 
quarks, we used the total energy density, including quarks, to find the initial tem- 
perature. The entropy increase associated with chemical equilibration of quarks 



was estimated to be about 15% at LHC in |jTT|. For SPS one has the problem of 
the long transit time. 



5 Conclusions 

We have in this article computed how the initial values of the energy and number 
densities of QCD matter produced in ultrarelativistic central heavy ion [A + A) colli- 
sions depend on A and ^/s. These values were converted to final observed numbers for 
dNch/drj and dEx/dr] by assuming adiabatic expansion. This assumption is justified 
by the observation that the initial values very closely satisfy the kinetic thermalisation 
condition e = 2.7nT for a dominantly gluonic system at all A, ^/s. 

The computation is technically one in perturbative QCD, leading order with the NLO 
corrections estimated by a constant i^-factor, and including nuclear shadowing, but 
also nonperturbative physics enters indirectly in that the computations are performed 
at one transverse momentum scale, the saturation momentum psat- This is determined 
by a geometric transverse saturation criterion and is supposed to effectively represent 
all scales, both those > psat but also those < psat- In time, the initial production is 
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assumed to effectively take place at r = 1/psat, instead of a gradual dumping of energy 
and entropy from the collision lasting until ~ 1/Aqcd- 

Even at the SPS, where one normally expects soft stringlike physics to dominate, 
this entirely hard effective approach gives rather reasonable numbers. Maybe this gives 
some credibility to the predictions for RHIC and LHC, where one expects the approach 
to work better. From a purely theoretical point of view one can think of corrections 
from many different sources. 
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